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ABSTRACT 



Aims. To improve the parameters of the HD 17156 system (peculiar due to the eccentric and long orbital period of its transiting 
planet) and constrain the presence of stellar companions. 

Methods. Photometric data were acquired for 4 transits, and high precision radial velocity measurements were simultaneously 
acquired with SARG@TNG for one transit. The template spectra of HD 17156 was used to derive effective temperature, gravity, and 
metallicity. A fit of the photometric and spectroscopic data was performed to measure the stellar and planetary radii, and the spin-orbit 
alignment. Planet orbital elements and ephemeris were derived from the fit. Near infrared adaptive optic images was acquired with 
ADOPT@TNG. 

Results. We have found that the star has a radius of R^ = 1 .44 ± 0.03Ro and the planet Kp = 1 .02 ± 0.08R j . The transit ephemeris is 
T, = 2 454 756.73 1 34 ± 0.00020 + • 21 .2 1 663 ± 0.00045 BID. The analysis of the Rossiter-Mclaughlin efi'ect shows that the system 
is spin orbit aligned with an angle /i = 4.8° ± 5.3° . The analysis of high resolution images has not revealed any stellar companion with 
projected separation between 150 and I 000 AU from HD I7I56. 

Key words, stars: individual: HD I7I56 ~ binaries: eclipsing - planetary systems - techniques: spectroscopic, photometric 



1. Introduction 

The discovery of transiting extrasolar planets (TESP) is of 
special relevance for the study and characterization of plane- 
tary systems. The combination of photometric and radial ve- 
locity measurements allows to measure directly the mass and 
radius of an exoplanet, and hence its density, which is the 
primary constraint on a planet's bulk composition. Dedicated 
follow-up observations of TESP during primary transit and sec- 
ondary eclipse at visible as well as infrared wavelengths al- 
low direct measurements of planetary emission and absorp- 
tion (e.g., ICharbonneau et al] (12007.) . and references therein). 
Transmission spectroscopy during primary eclipse in recent 
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years has been successful in characterizi ng the atmospheric 
chemi stry of several Hot Jup iters (e.g., ICharbonneau et all 
(I2002h and iTinetti et aL I (120071) ). Infrared measurements gath- 
ered at a variety of orbital phases, including secondary eclipse, 
have permitted the characterization of the longitudinal temper- 
ature profiles of nearby TESP. The quickly increasing amount 
of high-quality data obtained for TESP has provided the first 
crucial constraints on theoretical models describing the physi- 
cal structure and the atmospheres of gas and ice exoplanets. The 
detailed characterization of TESP ultimately is of special rele- 
vance to test several proposed formation and orbital evolution 
mechanisms of close-in planets. 

Th e planet HD 17156b, detected by iFischer et al.l 
(|2007|) (hereinafter F07) using the radial velocity method, 
was sh own to transit in front of its parent star by Barbieri et al.l 
(l2007h . Additional ph otometric measurements were presente d 
in follow-up papers by 'Gillon et al.' ('20081). iNarita et all (|2008|) . 
Irwin et al. (2008), Winn et al. (2008FThis planet is unique 
among the known transiting systems in that its period (21.2 
days) is more than 5 times longer than the average period for 
this sample, and it has t he larges t eccen tricit y (e = 0.67). 

ISchlesinieil (11910'). Ros siteii (1 1 924 and lMcLaughlinl (1 1 9241) 
showed that a transiting object, Uke a stellar companion, pro- 
duces a distortion in the stellar line profiles due to the partial 
eclipse of the rotating stellar surface during the event, and thus 
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an apparent anomaly in the measured radial velocity of the ob- 
served primary star In the case of close-in giant planets transit- 
ing solar-type stars, the amplitude of the Rossiter-McLaughlin 
effect ranges typically between a few and ~ 100 m s"', depend- 
ing on orbital period, stellar and planetary radii, and stellar rota- 
tional velocity. This is within the current instrument capabilities 
for the brightest transiting planets. Indeed, this effect has been 
previously observed on several TESP (e.g. Winn (2008) and ref- 
erences therein). The observation of the Rossiter-McLaughlin 
effect allows to measure the relative inclination angle between 
the sky projections of the orbital plane and the stellar spin axis. 
Another important result derived from these observations is the 
possibility to determine whether short orbital period transiting 
planets move in the same direction of the stellar spin, indi- 
cating the existence of strong dynamical interactions between 
the parent star and its planet. Measurements of the Rossiter- 
McLaughlin effect thus provide relevant fossil evidence of planet 
formation and migration proc esses, as well as dynamical interac- 
tions with perturbing bodies dMarzari & Weidenschillin j, |2002|) 
Up to now 8 of the 9 transiting planets for which the Rossiter- 
McLaughlin effect was measured are coplanar systems, as our 
Solar System (the orbital axes of all the planets and the Sun spin 
axis are aligned within a few degrees). This is compatible with a 
formation mechanism for close-in giant planets including migra- 
tion via tidal interactions with the protoplanetary disk. Only the 
XO-3 planet seems to be a non aligned system (Hebrard et al], 
120081) . 

The measurement of the Rossiter-McLaughlin effect is of 
special relevance for the HD 17156 system, given the very high 
eccentricity and orbital period, much longer than the other tran- 
siting planets. With a period of 21 days, HD 17156b is well 
outside the peak in the period distribution of close-in planets 
at about 3 days, possibly implying a different migration history 
with respect to the other transiting planets. An eccentricity as 
high as that of HD 17156b can be hardly explained by models 
of planet migration via tidal interactions with the protoplane- 
tary disk. Alternatively, high-eccentricity planets might be the 
outcome of a variety o f possible planet-planet dynamical inter- 
actions (Marz ari & We idenschilling, 2002). In these scenarios, 
high relative inclinations between the stellar rotation axis and the 
planet orbit plane are possible. An additional way to get large 
relative inclinations and eccentricities is represented by Kozai 
resonances with close stellar companion. Interestingly, a rather 
large planetary mass coupled with a short orbital period might 
be an indication for binarity, as t ypically high-mass short-pe riod 
pla nets occur in binary systems (:Desidera & Barbierill2007h . 

iNarita et"aLl (|2008) presented the first observations of the 
Rossiter-McLaughlin effect in HD 17156. They found an angle 
between the sky projections of the orbital axis and the st ellar ro- 
tation axis A -62 ±25°. However. ICochran et all (|2008|) did not 
confirm this claim, suggesting instead well-aligned axes. Such 
a discrepancy calls for additional high-precision RV monitoring 
during planetary transit. These are presented in this work, along 
with additional photometric observations and a critical revision 
of stellar and planetary parameters. To further constrain the ori- 
gin of the special properties of HD 17156b , we also searched 
for possible wide stellar companions using adaptive optics ob- 
servations. 

The overall organization of this paper begins in Sec. |2] with 
the description of high resolution spectroscopy data obtained 
with SARG@TNG. The following Sec|3]covers the analysis of 
the stellar parameters. Sec|4]presents the photometric data col- 
lected during several planetary transits. Next Sec|5]describes the 
analysis of the radial velocity and photometric data, and their re- 



Table 1. Differential radial velocities and bisector velocity span 
of HD 17156 obtained with SARG at TNG. 



BJD 


RV 


err 


span 


err 




m/s 


m/s 


m/s 


m/s 


2454398.58856 


-81.03 


5.76 






2454438.40203 


100.47 


4.55 


-25.47 


63.07 


2454438.41331 


82.40 


4.56 


-32.78 


62.90 


2454438.42460 


76.72 


4.55 


19.04 


60.32 


2454438.43608 


90.10 


4.43 


74.70 


59.30 


2454438.44737 


91.46 


3.98 


-5.58 


56.75 


2454438.45866 


76.28 


4.35 


64.69 


59.89 


2454438.47011 


65.86 


3.98 


47.68 


57.10 


2454438.48140 


60.08 


4.01 


40.49 


53.96 


2454438.49270 


30.74 


3.72 


119.82 


57.44 


2454438.50415 


32.97 


3.98 


216.98 


59.74 


2454438.51544 


20.63 


4.49 


108.08 


57.92 


2454438.52671 


18.78 


4.59 


74.54 


57.66 


2454438.53816 


4.29 


4.29 


82.66 


58.07 


2454438.54946 


-3.77 


4.13 


90.90 


57.79 


2454438.56075 


7.21 


4.93 


41.59 


63.36 


2454438.57218 


-3.84 


3.77 


58.41 


57.50 


2454438.58349 


-11.52 


4.53 


50.75 


60.76 


2454438.59478 


-31.71 


4.43 


45.47 


59.96 


2454438.62486 


-26.61 


4.74 


36.70 


64.19 


2454438.63615 


-43.01 


4.01 


18.82 


58.87 


2454438.64745 


-46.79 


4.23 


66.00 


65.44 


2454438.65894 


-54.19 


4.77 


-38.79 


65.29 


2454438.67023 


-55.35 


4.72 


3.19 


67.30 


2454438.68152 


-64.85 


5.32 


102.93 


70.23 


2454438.69306 


-71.94 


5.20 


73.98 


67.86 


2454438.70435 


-82.48 


5.08 


111.26 


72.47 


2454438.71584 


-88.29 


5.13 


-18.92 


72.31 


2454438.72713 


-92.59 


5.15 


94.27 


72.55 



suits. In Sec|6]we present the results for the search of additional 
stellar companions to HD 17156 . Finally in Sec|7]our conclu- 
sion are presented. 

2. High resolution spectroscopy 

We observed HD 17156 on 2007 December 3, including contin- 
uous monitoring (about 8 hours) during the transi t, with SARG, 
the high resolution spectrograph of the TNG dGratton et all 
l200lb . Observing conditions were not optimal, with seeing rang- 
ing from 1.3 "to 2.0 ". We obtained the stellar template (with- 
out the iodine cell) first and then started the uninterrupted se- 
ries of observations with the iodine cell. Exposure time of the 
spectra acquired with the iodine cell was 900 s, typically result- 
ing in S/N of about 80 per pixel. One additional spectrum was 
obtained on 25 Oct 2007. These spectra were reduced and ana- 
lyzed as those for the o ngoing planet search program with SARG 
(Desidera et alj, l2007h using the AUSTRAL code dEndl et all 
2000). Table[T]lists the radial velocities. 

The absolute radial velocity of HD 17156 was derived by 
cross correlating the stellar template acquired with SARG to 
a few suitable reference stars observed with the same set-up 
and with avai lable high-accuracy absolute radial velocity from 
,Nidever et al.k200Z) . It results in -3.15 + 0.20 km/s. 

3. Stellar parameters 

3. 1 . Spectroscopic analysis 

We have used the TNG/SARG template spectrum of 
HD 17156 to provide an independent assessment of its at- 
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mospheric parameters (Teg, logg, and [Fe/H]) with respect to 
the values reported by F07. 

Our methodology follows a standard procedure whose de- 
tails can be found in several works of the recent past (e.g. , 
Gonzalez & Lambert! [T996HGonzalez et al.Ll200ll:ISantos et all 



HD 17156 



1.0 



2004 . We briefly summarize it here. We initiall y selected a set 
of re latively weak Pel and Fell lines (see, e.g.. lSozzetti et al.l 
12004 , and references therein, for details on the line list), and 
measured equivalent widths (EWs) using the auto mated soft- 
ware A RES Q, made available to the community by ISousa et akl 
( I2007h . The EWs measured with ARES are then fed to the 2002 
version of the MOOG spectral synthesis code (Snedenl Il973h 
B, together with a g rid of Kurucz ATLAS plane-parallel stellar 
model atmospheres (lKuruczlll993h . 

The atmospheric parameters of HD 17156 are then de- 
rived under the assumption of local thermodynamic equi- 
librium, using a standard t echnique of Fe ionization bal- 
ance (see, e.g., I Santos et all l2004t ISozzetti et all 12004 and 
references therein). We obtained T^s = 6100 + 75 K, 
log^ = 4.1 ± 0.1, and [Fe/H]= -1-0.14 + 0.08, the for- 
mal errors on T^g and log g having been derived us i ng the 
pro cedure described in iNeuforge-Verheecke & MagainI (Il997h 
and iGonzalez & Vanturd (Il998l) . while the nominal uncertainty 
for [Fe/H] corresponds to the scatter obtained from the Fe I lines 
rather than the formal error of the mean. 

We also quantified the sensitivity of our iron abundance de- 
termination to variations of + lcr with respect to the nominal T^s, 
and log^ values, and found changes in [Fe/H] of at most 0.05 
dex, below the adopted uncertainty of 0.08 dex. 

With the aim of further testing the accuracy of the T^ff de- 
termination above, we have carried out a few additional consis- 
tency checks. For example, in Figure[T]we show the comparison 
of the observed H„ line profile in an archival Keck/HIRES spec- 
trum against four synthetic profiles for solar-metallicity dwarfs 
([Fe/H] = 0.0, \ogg - 4.5) from the Kurucz database. As is well- 
known, the Hq, line is very sensitive to changes in T^s, while 
relatively insensitive to changes in log; ? and [Fe/H] (see, e.g., 
ISantos et al.Ll2006tlSozzetti et al.Ll2007 l, and references therein), 
thus this exercise certainly helps to test the accuracy of the spec- 
troscopic Teff derived above. The results shown in Figure [1] in 
which a 10 A region centered on Hq. is displayed together with 
four calculated profiles for different Tf.g values, indicate rather 
good agreement with the estimate reported in Table |2] 



3.2. Age 



Based on isochrone fitting, F07 reported an age estimate for 
HD 17156 of 5.7^[ g Gyr, suggesting an old, slightly evolved 
F8/G0 primary. We performed an independent isochrone fitting 
using the set of isochrones of Girardi et al. (2000) and the soft- 
ware PARAM described in da Silva et al.l i2006B The input val- 
ues for PARAM are the parallax, the visual magnitude, [F e/H] 
and Jeff. We used the parallax from Ivan LeeuwenI (I2007h and 
for the metallicity and effective temperature we ran the code 
twice, once with the values from F07 and once with our esti- 
mate. The results of PARAM for the stellar age are 2.4 +1 and 
2.8 +1 Gyr, for the F07 and our parameters, respectively. These 
are only marginally compatible with the previous age estimate 
by F07. The stellar mass and radius are instead fully compatible 
(see below). 
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Fig. 1. The portion of the spectrum of HD 17156 around H^ and 
synthetic profiles for four different T^g. A T^ff close to 6000 K is 
suggested by the fit of the wings of the line. 



Other indirect age indicators confirm an age of a few Gyr. 
The low level of Ca II H&K chromospheric activity suggests 
an age of about 6 Gyr (Fischer et al .l | 20p7l). The lack of X- 
ray emission from ROSAT (Vogesetal., 2000) yields an up- 
per limit of logLx < 28.7. This in turn implies an age older 
than 1.6 Gyr, using the age- X ray emission calibration by 
iMamaiek & Hillenbrandl (l2008l) . 

To obtain an additional age estimate and to investigate pos- 
sible chemical peculiarities of HD 17156 with respect to other 
planet hosts with similar physical properties, we have measured 
its lithium (Li) abundance. 

Figure |2] shows a spectral synthesis of a 10 A region cen- 
tered on the Li A - 6707.8 A line in an archival Keck/HIRES 
spectrum of HD 17156 , and using the atmospheri c parameters 
derived from the Fe-line analysis and the line list of Reddv et"aLl 
(l2002h . In the figure, the observed spectrum is compared to 
three synthetic spectra, each differing only in the assumed Li 
abundance. We find a best-fit value of log e(Li) a; 2.80 for 
HD 17156 . We then infer a rather old age for the star of f > 2 
Gyr, based on the average Li abundance curves as a function 
of effective temperatu re for clusters of different ages reported 
bv lSestito & Randichl (f2005). 

The measured Li abundance for HD 17156 does not ap- 
pear peculiar when compared to that of sub- sample s of nearby 
planet hosts with similar T^g (llsraelian et al.L l2004t iGonzalea, 
2008). To further investigate these issues we will present in a fu- 
ture paper a more detailed study of the elemental abundances in 
HD 17156. 



3.3. Stellar mass and radius 



http : //verdi . as . utexas . edu/moog . html | 
http : // stev . oapd . inaf . it/param| 



F07 provided mass and radius estimate from isochrone interpo- 
lations (mass 1.2 ±0.1 M©, radius 1.47^" }^ R©). Our isochrone 
fit (see Sect. O gives 1.21 +0.03 and 1.19 +0.03 Mq when 
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Fig. 2. Portion of HD 17156 spectrum centered on Lithium 
6707.8 A line. Overplotted the results of the spectral synthesis 
for three different lithium abundance. 



assuming Teff and [Fe/H] from the F07 and our analysis respec- 
tively, and same radius 1 .35 +0.08 R© 

The comparison of the results of the two fits show that they 
are compatible within error bars, but the best value are slightly 
different. The origin of these differences is not univocal, and to 
understand this problematic independent measurements of the 
star mass and radii are needed. 



3.3.1. Monte Carlo experiment. Method 

For this purpose we adopted an approach based on a Monte Carlo 
experiment. In each realization we generated a set of observed 
stellar properties. Using using calibrations from the literature we 
obtained the corresponding values for mass and radius. From the 
resulting distributions we obtained the most probable values and 
relative errors for the mass and radius. 

In more detail, we created 10'' different synthetic systems 
where we generated Gaussian distributions for the parallax, the 
V and K magnitudes, and T^ff, using as standard deviations their 
respectively error bars. The input parameters and relative their 
standard deviations are reported in Tab. |2] In the same way, we 
generated values for the bolometric correction and the bolomet- 
ric magnitude of the Sun. For each synthetic set obtained this 
way we calculated the absolute magnitudes, luminosity, radius, 
mass, density and gravity. 

The V magnitude was obtained from Simbad, K magni- 
tude from 2MASS a fter conversion to the Bessel-Brett sys- 
tem dCarpenteii 1200 ll) . and the T^ff from our spectroscopic de- 
termination. For t he parallax we adopt ed the recently revised 
Hipparcos value jvan Leeuwenl l2007h . which indicates that 
HD 17156 lies at 75 pc from the Sun, 3 pc closer than the previ- 
ous estimate. The bo lometric correction was set to B.C.= -0.03 
dOirardi et all l2000l) while for the solar bolometric magnitude 
we used the value M^o/ o = 4. 77 ±0.02. 

Th e stellar radius w as obtained using the Teff-K mag calibra- 
tion of iKervella et alj (|2004|) . and using the Stefan-Boltzmann 
law. The stellar mass was calculated from various mass- 
luminosity relations ( MLR) namely: i) the classical MLR L cc 



M^^; ii) the MLR of iMalkovl mO% using absolute magni 



tude and stellar luminosity; iv) the MLR of Henry (2004). 
Density and gravity were estimated directly from the radius and 
mass assuming a spherically symmetric star 



3.3.2. Monte Carlo experiment. Results 

The distribution of absolute magnitudes in V band peaks at 3.7 
(Fig. [3] upper left), while the luminosity is peaked at 2.5 Lq 
(Fig. [3] upper rig ht). These valu es are slightly different from the 
ones obtained by iFischer et all and the source of these differ- 
ences is ascribed only to the difference in the adopted parallax. 

The resulting distributions for the radius are shown in Fig. [3] 
(middle left). The two relations used provide similar results 
(Rs ~ 1 .4R0) for the b est value and also the shape of the distri- 
butions is very similar. [Fischer et al] suggest a slightly larger ra- 
dius. Also in this case, the difference originates from the change 
in the adopted parallax. 

In Fig. [3] (middle right) we present the mass distributions ob- 
tained with the MLRs. Using the relation of Henry. (2004) we ob- 
tain the highest mass (M ~ I.28M0). However, we note that this 
relation was originally derived from parameters of close binary 
stars. Malkov (2003) demonstrated that this kind of relations do 
not well describe single stars. In order to avoid this problem 
we adopted the MLRs of lMalkovl (|2007[) obtained on detached 
main-sequence double-Uned eclipsing binaries. These relations 
are also valid for slightly evolved stars like HD 17156 (almost 
all the stars used for deriving these relations are also slightly 
evolved. O. Malkov, private communication). We obtain best 
values for the mass between 1.2 and 1.24 M© , the first obtained 
using MLR(Mv) and the second using MLR(L). The classical 
MLR (M cc L**^) provides M -1.22 Mq. These results are con- 
sistent with the values estimated by F07 (1.2 + 0.1 Mq). 

Finally, the resulting distributions for the gravity and the den- 
sity are portrayed in the lower panels of Fig. [3] The mean values 
are \ogg - 4.22 andp = 0.58 g/cnr'. 

The results of this experiment show a fairly good agreement 
with the estimate of Rs based on isochrone fitting (both our 
and the one by F07). We conclude that the two independent ap- 
proaches based on isochrone fitting and the use of scaling rela- 
tions provide consistent results. 

In the following we do not adopt a value for the radius, be- 
cause we want to determine independently its value from the 
light-curve fits. Moreover, we fix the value of the mass to the 
value of the weighted mean of our mass estimation (not using 
the Henry MLR results) M = 1.24+0.03 Mq. We summarize in 
Table|2]all the data relative to this Monte Carlo experiment. 



3.4. Rotational velocity 

F07 derived vsin/ = 2.6 + 0.5 km/s for HD 17156 . Our tem- 
plate spectrum is suitable for an independent measurement of 
this quantity. We discuss here three different methods adopted 
for the measurements of vsin /, based on our template spectra. 

With the first method we derive vsin / by the Fast Fourier 
Transform analysis of the star's absorption profile (see Fig|4|i. 
To determine the vsin/, the observed profile of a stellar ab- 
sorption line is made symmetric by mirroring one of its halves, 
with the purpose to reduce the noise of the FFT. A new pro- 
file is calculated by the convolution of a macroturbulence pro- 
file (Gaussian) and a rotational one, to compare the FFTs of 
the symmetric and the calculated (model) profiles. The vsin/ 
value of the rotation profile, is set as variable parameter until 
the first minimum of the FFT from the calculated profile coin- 
cides with the minimum of the FFT from the symmetric one. 
The value of vsin / for HD 17156 was determined considering 
poss ible values of macroturbu lent velocity iv,„ac) from B-V and 
Teff dValenti & Fischeii l2005l) . and we obtained a vsin / ranging 
from 1.8 to 2.8 km /s. 



M. Barbieri et al.: Characterization of the HD 17156 planetary system 



5 



Table 2. Upper panel Input parameters for the Monte Carlo 
experiment and best values. Lower panel Kinematical properties 
and galactic orbit parameters. 



Monte Carlo experiment 
Input parameters 



parallax 


13.33± 0.72 


mas 


mag V 


8.172 ±0.031 


mas? 


mag K 


6.807 ± 0.024 


mag 




6100 ± 75 


K 


B.C. 


-0.03 ± 0.02 


mag 


Output parameters 




Mboi,o 


3.69 ±0.12 


mag 


Mv 


3.73 ±0.12 


mag 


L 


2.68 ±0.28 


Lg 


R (Stefan-Boltzmann) 1.49 ±0.09 


Rq 


R (Kervella) 


1.45 ±0.07 


Rq 


M (Malkov,Mv) 


1.21 ±0.04 


Mo 


M (Malkov,L) 


1.25 ±0.04 


Mo 


M (M oc L'*-^) 


1.22 ±0.04 


Mo 


M (Henry,Mv) 


1.29 ±0.03 


Mo 


log g (mean) 


4.21 ±0.05 


cgs 


p (mean) 


0.57 ±0.10 


g/cm^ 


Kinematical properties 




RV 


-3.15 ±0.2 


km /s 




91.14 ± 0.49 


mas/yr 




-33.14 ±0.56 


mas/yr 


U 


0.6 ± 0.2 
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Fig. 3. Probability density function of the stellar parameters ob- 
tained with the Monte Carlo experiment. Upper left: Absolute 
magnitude Upper right: Luminosity Middle left: Radius: thick 
line results using Stefan-Boltzmann law, thin line results using 
the calibration of Kervella(reff,MK) Middle right: Mass, from 
left to right results using: Malkov MLR(Mv), M oc Lf'-^, Malkov 
MLR(L), and Henry MLR(Mv) Lower left: log g Lower right: 
Density 



The second method that we used consists in obtaining the ro- 
tational velocity by means of a suitable calibration of the FWHM 
of the cross-correlation function against the B-V color This re- 
lation was derived for all the stars in the SARG planet search 
survey, and it was calibrated into vsin / using stars with known 
rotational velocity from the literature. Using the B-V from the 
Tycho catalog converted to the Johnson system (B-V=0.632), 
the resulting vsin / = 3.2 + 1 km /s . 

Finally, using MOOG we synthesized a number of isolated 
Pel lines in the template spectrum. From these we measured 
vsin / = 3.0 + 0.5 km /s. 

The values obtained with the three methods suggest a range 
of values for vsin / compatible with the measurement of F07. The 
measurement of vsin / from the analysis of Rossiter-Mclaughlin 
effect will be presented in Section[ 



HD17156 




Log Fourier cycles (s/km) vsini = 2.83000 (kr^/s) > 



= 4.4553B (km/s) 



3.5. Galactic orbit 

The measurements of the absolute radial velocity given in 
Section |2] toget her with the revised parallax and proper motion 
from Hipparcos'v an LeeuwenI (l2007h . allow one to calculate the 
space velocity with respect to the Local Standard of Rest and 
the galactic orbit of HD 17156 . Space v elocities are calculated 
follow ing t he procedure d elineated by iJohnson & SoderblomI 
(Il987l) and iMurravl (1 19891). adopt i ng the value of standard so- 
lar motion of iDehnen&BinnevI (Il998h (with U positive to- 
ward the galactic anticenter). The calculations yield {U, V, W) - 
(h-0.6, -1-26.1, -22.8) km /s. The galactic orbit of the star is ob- 
tained integrating the equations of motion of a massless parti- 



Fig. 4. Estimation of vsin ; through the FFT Three profiles were 
analyzed: the observed profile (solid line), the symmetric profile 
(dashed line) mirroring one half of the profile, and a theoretical 
profile (dotted line). The FFT corr esponds to a case conside ring 
Teff = 6079 K and Eq. 1 given by (IValenti & Fischeiil2005h . 



cle in the potential described by lAllen & SantillanI (Il99lh . The 
equations of m otion are solved using the RADAU integrator 
('Everhart'.'1985) assuming that the rectangular galactocentric co- 
ordinates of the Sun are (X^, Yq,Zq) = (8.0, 0.0, 0.015) kpc, and 
that the local circular velocity is 220 km s ' . We compute 1 000 
orbits each time varying randomly the initial coordinates and ve- 
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locity of the star within the error bars and integrating for a time 
of 4 Gyr (^15 full orbits around the Galactic Center). The UVW 
spatial velocity and the mean values of the computed orbits are 
reported in Table |2] The mean galactic radius provides an esti- 
mate of the galactocentric distance of the star at the moment of 
its formation. The value for HD 17156 , Rmed = 9.5 kpc, implies 
that the star spends most of its time in regions outside of the solar 
circle. 

4. Photometric observations 

Photometric observations of the transit of 2007 December 3 were 
obtained with various telescopes. On 2007 December 2 we used 
three medium-class telescopes (Asiago 1.82 m and OHP 1.20 
m) and a number of small telescopes, including six 30-40 cm 
amateur-operated telescope all located in continental Europe as 
well as the Telast 0.3-m telescope in the Canary Islands. Weather 
conditions across continental Europe were not optimal. 

Observations at Asiago started under photometric condi- 
tions, transit ingress was observed but observations ended at 
December 04 UT 02:30, due to the presence of clouds, which 
prevented the observations of the third and fourth contacts. OHP 
observations were performed under variable sky conditions due 
to intermittent clouds and veils, good observing conditions were 
achieved only during the transit time window, with only a small 
coverage of the Out Of Transit (OOT) flat part of the lightcurve 
before the first contact. Observations with Telast were obtained 
under normal sky conditions and were performed throughout the 
night. Amateur observations were carried out by six observato- 
ries spread over central and northern Italy. Observing conditions 
suffered from clouds and veils as the others European sites in- 
volved in the campaign, and the full transit was successfully ob- 
served by four telescopes, the remaining two observatories ob- 
taining data only for the ingress phase of the transit. 

Three additional attempts to coordinate transit observations 
of HD 17156b were carried out on 25 December 2007 and on 
September 25 and October 17 2008. Unfortunately, due to bad 
weather conditions, there was only one useful observation for 
each date. 

All observations were obtained in R band except the two ob- 
tained by B. Gary that was acquired in white light; the charac- 
teristics of the different telescopes are summarized in Table |3] 

4.1. Data reduction 

The Asiago, OHP, and Telast raw images were calibrated using 
flat field and bias frames. The resulting images were analyzed 
with IDL routines to perform standard aperture photometry. The 
center of the aperture was calculated using a Gaussian fit, and 
the aperture radius was held fixed for each set (in the range 15- 
20 pixels). The sky background contribution was removed after 
an estimation The brightest non-variable stars in the field were 
measured in the same way, and a reference light curve was con- 
structed by adding the flux of these stars. The target flux was 
divided by this reference to get the final normalized light-curve. 

The images obtained by amateur telescopes were reduced 
bias subtracted, and dark flat-field corrected using commer- 
cial software. Aperture photometry of HD 17156 was then per- 
formed using IRIS Q and adopting a fixed aperture equal to 2 
times the stellar FWHM. The sum of the flux of the brightest 
stars in the field was used as reference for building the normal- 
ized light-curve. 

* ihttp : //www ■ astrosurf ■ com/buil/iris| 



Table 3. Summary of photometric observations of 
HD 17156 performed on 2007/08. The last column report 
the midtransit times for each date obtained from the fit in 
Sec 1121 



site/Observer 


diatnctcr 


rms OOT 


date 


T, 




m 


mtnag 




BJD 


Telast 


0.30 


4.2 


2007 September 10 


2454353.61300 ± 0.0200 


Gasparri 
Lopresti 


0.20 
0.18 


5.3 
7.0 






Asiago 


1.82 


2.8 


2007 December 03 


2454438.47450 ± 0.0005 


OHp" 


1.2 


2.0 






Telast 


0.3 


11.2 






Obs. Cavezzo 


0.4 


6.1 






Lopresti 


0.18 


10.5 






Obs. Univ.Siena 


0.3 


8.4 






Obs. IVIt. Baldo 


0.4 


7.3 






Papini 


0.3 


6.9 






Vallerani 


0.25 


5.9 






Gary 


0.3 


5.6 


2007 December 25 


2454459.69087 ± 0.0032 


Gregorio 


0.3 


3.8 


2008 September 25 


2454735.51351 ±0.0027 


Gary 


0.3 


8.9 


2008 October 17 


2454756.73134 ± 0.0024 



The final light-curve for each telescope was corrected for dif- 
ferential airmass and residual systematic effects dividing them 
by a linear function of time to the region outside the transit. The 
photometric error on each point of a lightcurve was calculated 
as the rms over an interval of 30 minutes (the timescale of the 
ingress/egress phase). The typical rms of the OOT lightcurves 
are reported in Table |3] while the complete photometric dataset 
will be available in electronic format at CDS. 

The whole dataset consists of ~ 7 000 photometric points. 
We used these data to perform a global analysis of the planetary 
transit. 

For the light curve fitting we used all the lightcurves that 
we have collected without performing data binning. In Fig |5] are 
portrayed the light curves used in this study, folded with the best 
orbital period from the fit. For displaying purposes the combined 
light-curve of the 15 lightcurves is shown in Fig.|6] This com- 
bined light curve was obtained using a bin width of 90 s, the 
OOT has an rms of 0.0016. 



5. System parameters 

We performed the analysis of the HD 17156 system in three 
steps. First, using the radial velocities presented in Table 1 along 
with other published RV values : F0 7 (2 datasets: Keck+Subaru), 
[Narita et al. {2008), ICochran et all (|2008) (2 datasets: HET -i- 
HJST), we have derived a new spectroscopic orbital solution for 
HD 17156 . A full Keplerian orbit of five parameters: the radial 
velocity semi-amplitude K^y, the time of periastron passage Tp, 
the orbital period P, the orbital eccentricity e, and the argument 
of periastron o) was adjusted to the data. Second, we have car- 
ried o ut a fit to all the ligh- curve that we have obtained and also 
to the iBarbieri et alj (|2007|) datasets using the e, and the oj ob- 
tained from the orbital solution. In this scheme, the adjustable 
parameters are the ratio of the radii k - Rp/Rs their relative 
sum {Rs + Rp)/a, the orbital inclination /, the midtransit time Tc, 
and the orbital period P. The values derived from the light-curve 
analysis were then used to determine, through the analysis of the 
Rossiter-McLaughlin, new values of vsin / and of the angle A be- 
tween the equatorial plane of the star and the orbital plane of the 
planet. 

The modeling of the transit lightcurve and Rossiter- 
McLaughlin effect was carried out using the analytical formu- 
lae provided by Gim enez (2006a) and Gimenez (2006b). The 
mathematical basis for the description of the two effects is the 
same, i.e. the lKopall(ll977h theory of eclipsing binary stars. This 
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Fig. 6. Combined light curve of all photometric data folded with the orbital period, along with the best fitting model. The time 
interval between each point is 90 s. 



fact warranties an internal consistent description of the observed 
data. 



5.1. Orbital radial velocity analysis 

In order to derive the stellar spectroscopic orbit using the com- 
bined set of radial velocities mentioned above we used only the 
OOT measurements in all datasets. Observations in the night of 
the transit are valuable to this goal because of the steep RV slope 
(about 23 m/s/hr). We use a downhill simplex algorithm to per- 
form the RV fit to the six datasets, including the zero point shifts 
between the datasets as free parameters. A stellar jitter of 3 m/s 
was added in quadrature to the observational errors F07. 

The best-fit solution has a value of reduced - 1.13, and 
the results are in close agre ement with the dis covery paper F07 
and its subsequent analysis dlrwin et aUl2008h . Uncertainties in 
the best fit parameters were obtained exploring the grid with 
an adequate resolution. The orbital solution and relative param- 
eter uncertainties are presented in Table |4] In Figure |7] we show 
the phased radial velocity curve with the best-fit model. Using 
the value of primary mass provided in Section 13.31 and its un- 



certainty, the resulting minimum mass for the planet is m sin / = 
3.21 + 0.08 Mj, and the semi-major axis is a = 0.1614 + 0.0010 
AU. 



5.2. Photometric analysis 

We used the description of lGimene3 ( l2006al) to analyze the light- 
curves obtained in Sect. 14.11 In our model we allowed to vary 
the ratio of the radii, the phase of first contact, the time of tran- 
sit center and the orbital inclination. We fixed the limb darken- 
ing coefficients to the values corresponding to a star with simi- 
lar temperature and metallicity to HD 17156 from Claret (200^ 
tables. For the R band the adopted limb darkening coefficients 
were : u+ = Ua + ui, = 0.6323 and m_ - - Ub = -0.0655. 
The eccentricity, and the longitude of periastron were held fixed 
to the best fit values obtained from the RV analysis (Sec. 15.1b . 
Errors were estimate d using the bootstrap scheme described in 
(lAlonso etal.L 120081) 

The results of of the analysis of the 15 datasets are collected 
in Table |4] Using the third Kepler's law we obtain for the stellar 
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Table 4. Parameters of the HD 17156 system. 
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Fig. 5. Mosaic of the differential light curves obtained dur- 
ing transits of HD 17156b . obtained with several tele- 
scopes. In each box the horizontal axis is the photometric 
phase and the vertical axis is the relative flux, along with 
the best fit model. From left to right column and from top 
to bottom: Tla =Almenara, Tl b^Gasparri and Tlc=Lopresti 
datasets from flie lBarbieri et all |2007). T2a=OHP, T2b=Telast, 
T2c=Castellani, T2d=Asiago, T2e=Lopresti, T2f=Marchini, 
T2g=Nicolini, T2h=Papini, T2i=Vallerani. T2a to T2i light- 
curves were collected on 2007 December 03. T3a=Gary 2007 
December 25, T4a=Gregorio 2008 September 25, T5a=Gary 
2008 October 2008, 



and planetary radii Rs = 1.44 ± 0.08 Rq and Rp = 1.02 + 0.08 
Rj. 

The histogram of the residuals of the light curves (Fig.© has 
a gaussian shape with standard deviation of 0.0062. 

The result s are consistent wi t h pre v ious determinations 
dBarbieri et all 120071: llrwin et all 120081: iNarita et all l2008t 
iGiUon et all 120081) . nevertheless the values of incUnation and 
stellar radius show larger deviations with respect to values pre- 
sented by other authors^ 

In their analysis [irwin et alj (l2008h and lNarita et alj (|2008|) 
fixed the st e llar r adius to the value proposed by F07, while 
iGiUon et alJ (12008!) obtained directly the radius from their anal- 
ysis. The origin of the discrepancy on the stellar radius might 
lie in the different values for the inclination, because its value 
controls the transit duration and the relative sum of the radii. For 
confirmation we repeated the fit with only the published light- 





Orbital parameters 




p 
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Fig. 7. Upper panel: radial velocities of HD 17156 phased to the 
best-fit orbital solution. Bottom panel: residuals from the orbital 
solution. 



curve of iGillon et al.l (l2008h and keeping the limb-darkening 
coeflicients for the B band fixed to 0.7989 
and - Ua - Ub - 0.3019. The results are the following: 
k = 0.0729 + 0.0031, 01 = -0.00316 + 0.00023, / = 87.9 + 0.1, 
T, = 2454438.48372 + 0.00053 BJD, Rs = 1.44 + 0.07 R^ and 
Rp = 1.02 + 0.07 Rj. These results are very close to the results 
of our previous fit. 
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Fig. 8. Histogram of the lightcurve residuals. Overplotted the 
best Gaussian fit. 



We note that a transit model that does not take into account 
the non-zero eccentricity might lead to erroneous results in the 
orbital inclination a nd thus also the ste llar radius (see, for in- 
stance, section 3.2 in lAlonso et al.ll2008h . 

5.3. Rossiter-McLaughlin effect analysis 

The analysis of the TNG RV data obtained du ring the transit wa s 
performed using the formalism developed bv lGimene3 ( l2006"bl) . 
We allowed to vary vsin / and A and we fixed the values of K, 
P, e, cl>, k, i, Tc, 6i to the best values obtained from RV and 
photometry analysis and reported in Table |5] Fig|9]presents the 
best fitted model to the data. 

The best fitted values to the RV orbital residuals of SARG 
are vsin; = 1.5 + 0.7 km /s and A = 4.8° + 5.3°. The value of 
vsin i agrees with the values determined by F07 and by our anal- 
ysis of the stellar spectra. A is consistent with zero, indicating 
that the eclicptic plane of the planet is closely aligned with the 
equatorial p lane of the star. Th is value of A does not confirm 
the claim o f iNarita et alj (l2008l) for a large misalignment in this 
system, but rather ag rees with the relative alignment obtained by 
ICochran et al.l(l2008l) . Moreover, the two groups find very differ- 
ent values for vsin To study the nature of this discrepancy we 
repeated the fit on their datasets independently, and the results 
of the fits are summarized in Table|5] The results obtained using 
the HET dataset, in spite of their good precision, do not pro- 
vide strong constraints due to their partial coverage of the transit 
and the fact that the zero point in the OOT data cannot be esti- 
mated correctly. Instead, the fit to the HJST data are in excellent 
agreement with the determination obtained with SARG. Finally, 
the Narita et al. dataset provides a vsin / in good agreement with 
previous determinations, and a value of A that formally points 
toward the occurrence of some misalignment. These results in- 
dicate also that our adopted description is consistent with the one 
used by Narita et al. however, the intrinsically lower precision of 
their RV data makes these results not significant. 

We also measured line bisectors and bisector velocity span 
using the technique developed by iMartmez Fiorenzano et alj 
(l2005h and looked fo r changes in the line profiles caused by the 
planetary transit (see'Loeill et et al.ll20 08'). We do not detect sig- 
nificant variations (Fig.fTOli. This is not unexpected considering 
the typical signal to noise ratio of our spectra and the low ampli- 
tude of the Rossiter signature. 




I i 1 I CO I I i it! i 
9 i ■ i i ! i CSi 
] M I i M lii » 



0.005 
Photometric ptiase 



60 



50 
40 
30 - 
20 
10 



TNG 

OAO ^ 

HJST ( -■- 

HET : □ 











■ 


-■ 


■ 


V 








■ 


1* 




• 






■ ■ 












9 ( 





















0.005 
Photometric phase 



Fig. 9. Upper panel: residuals of radial velocities of 
HD 17156 phased to the best-fit orbital solution with the best 
fitted Rossiter-Mclaughlin effect overlayed. Bottom panel: resid- 
uals from the Rossiter-Mclaughlin effect. 

Table 5. Results of the Rossiter-McLaughlin modeling of all 
datasets. (OAO: Narita et al dataset). 



telescope 


date 


vsin i 


A 


x' 






km /s 


deg 




OAO 


17/11/2007 


1.8+ 1.5 


-22.6 ± 20.3 


0.7 


TNG 


3/12/2007 


1.5 ±0.7 


4.8 ±5.3 


1.4 


HJST 


25/12/2007 


2.2 ± 1.0 


0.8 ± 6.4 


0.9 


HET 


25/12/2007 


1.0 ±5.0 


30.1 ± 25.6 


0.3 


TNG+HJST 




1.6 ± 1.0 


3.9 ± 5.4 


1.2 



6. Clues on additional companions 

HD 17156 was obs erved with A dOpt@TNG, the adaptive op- 
tics module of TNG (ICecconi et al.. 2006) . The instrument feeds 
the HgCdTe Hawaii 1024x1024 detector of NICS, the near in- 
frared camera and spectrograph of TNG, providing a field of 
view of about 44 x 44 arcsec, with a pixel scale of 0.0437 
"/pixel. Plate scale and absolute detector orientation were de- 
rived in a companion program of follow-up of binary systems 
with lo ng term radial veloc ity trends from the SARG planet 
search jDesideraet al.Ll2007h . 

Series of 15 sec images on HD 17156 were acquired on 3, 
18 and 23 October 2007 in Bry intermediate -band filter Images 
were taken moving the target in diff'erent positions on the de- 
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Fig. 10. Bisector velocity span (BVS) from line bisectors of 
HD 17156 spectra. Upper panel displays BVS vs. RV. Lower 
panel shows the line bisectors from all spectra, where the hori- 
zontal lines enclose the top and bottom zones considered to com- 
pute the BVS. 



lector, to allow sky subtraction without the need for additional 
observations, and each night at three different field orientations 
to make it easier to disentangle true companions from image arti- 
facts. The target itself was used as reference star for the adaptive 
optics. Observing conditions were poor on the night of October 
3, and rather good on the nights of 18 and 23 October, when we 
obtained a typical Strehl Ratio of about 0.3. 

Data reduction was performed by first correcting for detector 
cross-talk using dedicated routines [j and then performing stan- 
dard image preprocessing (flat fielding, bad pixels and sky back- 
ground corrections) in the IRAF environment. Individual images 
taken at a fixed orientation were shifted and coadded. 

The successive analysis was optimized for the detection of 
companions in different separation ranges. At small separations 
(from about 0. 15 to 2 arcsec) we selected the two best combined 
images taken at different field orientations on 2007 Oct 18. They 
are shown in Fig. [TT] These two sets of images are character- 
ized by similar patterns of optical aberrations, and therefore, 
considering their difference, most of the patterns cancel out in 
difference images (Fig. [TT]), improving significantly the detec- 




Fig. 11. Images of HD 17156 . Upper panel: image of 
HD 17156 obtained with ADOPT@TNG. A similar quasi-static 
speckle pattern can be seen. Lower panel: difference between 
two images taken with two different field orientations, which al- 
lows to improve significantly detection limits in the inner re- 
gions. The field of view shown is 4.3x4.3". The images are dis- 
played in linear gray scale. 



http://www.tng.iac.eS/instruments/nics/files/crt_nics7.f 



tion limits (angular differential imaging, ' Marois et al.l (|2006')). 
In the differential image, a true companion is expected to show 
two peaks, one positive and one negative, at the same projected 
separation from the central star and position angle displaced by 
20° (the rotation angle between the two sets of images in our 
case). For detection at separations larger than about 2 arcsec, we 
summed all the images after an appropriate rotation, obtaining a 
deep image over a field of about 10x10 arcsec. 

No companion was seen in both the differential image at 
small separation and in the deep combined image within 10 arc- 
sec. The limit for detection was fixed at peak intensities 5 times 
larger than the dispersion over annuli at different radial separa- 
tion. The results, both for the differential image and the deep 
composite image are shown in Fig. [12] 

The contrast limits derived above were transformed into lim- 
its on companion masses using the mass-luminosity relation 
by iDelfosse et alj (l200d ). and projected separation in arcsec to 
AU using the Hipparcos distance to the star (Fig. [T3] ). A main- 
sequence companion can be excluded at a projected separation 
between about 150 and 1000 AU (the limit of image size). At 
such separations only brown-dwarf or white-dwarf companions 
are compatible with our detection limits. At smaller separation, 
detectability worsens quickly, and only stars with mass larger 
than about 0.4 M© can be excluded at projected separations 
closer than ~ 50 AU. 
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Fig. 12. Detectability limits for companions of difference mag- 
nitude around HD 17156 as a function of the projected separa- 
tion in arcsec Continuous line: limits on the difference image. 
Dotted line: limits on the composite deep image. 
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Fig. 13. Detectability limits for stellar companions around 
HD 17156 as a function of the projected separation in AU. 
Continuous line: limits on the difference image dotted line: lim- 
its on the composite deep image. 



The residuals from the radial velocity orbital solution do not 
suggest the occurrence of long term trends. This places further 
constraints on the binarity of the target. However, the timespan 
is rather short, and the continuation of the radial velocity moni- 
toring is mandatory for a more complete view. 

The available astrometric data from Hipparcos does not show 
evidence for stellar companions either (no astrometric accelera- 
tion within the timespan of the Hipparcos observations and no 
significant differences between 



7. Summary and Discussion 

In this paper, we studied the characteristics of HD 17156 and its 
transiting planets. Stellar parameters (mass, radius, metallicity) 
agree quite well with the previous study by F07. 

Our measurement of the stellar radius of HD 17156 obtained 
through the analysis of the transit light-curve is the same as the 
one obtained usin g the Stefan- Boltzmann law or the Kervella 
calibration (Tabl2ll lGillon et all (12008 ) obtained a radius of 1.63 



±0.2 Rq, which is only marginally compatible with our estimate 
1.44 ±0.08 Rq. On the one hand, to explain such a large stellar 
radius and the observed visual magnitude, it would be necessary 
to add ~ 0.3 mag of interstellar absorption, that at the distance 
of HD 17156 is not realistic (suggesting a mean extinction of 4 
mag/kpc) because HD 17156 is located well inside of the Local 
Bubble, where no strong absorption is present, and the maximum 
expected absorption is few hundreds of m ag. On the other hand, 
also the comparison of the iGillon et"ai1 (j2008) radius estimate 
with stellar models does not appear satisfactory: it is not pos- 
sible to find model with a radius that agrees with the observed 
temperature and metallicity of HD 17156 . We conclude that the 
determ ination of the stell ar radius, and by inference planetary 
radius. lOillon et al] ( l2008l) is overestimated by 15%. 

For a planet of 3Mj and an age of ~ 2 Gyr, theoretical models 
of planet evolution (Baraffe et al., 2008) predict a radius ranging 
between 0.9 and 1.1 Rj as a function of chemical composition 
of the planet. Our determination of the radius of HD 17156b is 
Rp - 1.02 ± 0.08 Rj. This is in excellent agreement theoretical 
expectations. Thus, the strong tidal heating effects on the planet 
do not ap pear to c ontribute to significantly inflate its radius. 

iFortnev et all (l2007i) suggested that HD 17156b , due to its 
large orbital eccentricity, can change its spectral type during the 
orbit from a warm pL type at apoastron to a hotter pM type at pe- 
riastron. Their models suggest that for pM class planets the ob- 
served radius at R band could be 5% larger than the radius mea- 
sured at B band (due to the increased opacity of TiO and VO in 
the R band). Comparing our radius measurements i n the R band 
and the measurements obtained in the B band by iGillon et al.l 
(I2OO8), we find identical results. This result is however not sig- 
nificant because of the large error involved in radius measure- 
ment. In order to obtain a significant difference of the radius in 
B and R band each measurement should be more accurate than 
0.03 Rj. 

Our RV monitoring of HD 17156 during the 2007 December 
3 transit does not confirm the misalignme nt between the stel - 
lar spin and planet orbit axes claimed by iNarita et"aLl (|2008|). 
but it a grees instead with the opposite finding by Cochran et al.l 
(l2008l) . We think our results are more robust because the other 
most accurate dataset (HET in iCochran et all l2008h does not 
cover the full transit, leaving some uncertainties in their Rossiter 
modeling. We then conclude that the projection on the sky of the 
stellar spin and planet orbit axes are aligned to better than 10 
deg. 

Therefore HD 17156 joins most of the other exoplanet sys- 
tems with available measurements of the Rossiter-McLaughlin 
effect in being compatible with coplanarity. The only possi- 
ble exception is rep resented by the XO-3 system, for which 
Hebrard et al. I (I2008h found indications for a large departure 
from coplanarity (~ 70°). However, as acknowledged by the au- 
thors, this result should be taken as preliminary, because of the 
possibility of unrecognized systematic errors for observations 
taken at large airmass and with significant moonlight contami- 
nation. 

Our result confirms that large deviations from coplanarity 
between stellar spin and planet orbit axes are at most rather 
rare. Such a rarity had already been established at a high con- 
fidence level for the "classical" Hot Jupiters in short-period cir- 
cular orbits. For massive eccentric planets the situation is less 
clear: HD 147506b (M = 8.6 Mj , P = 5.6 days, e = 0.52) and 
HD 17156 (M = 3.2 Mj , P = 21 days, e = 0.67) have projected 
inclinations below 10° while the possible detection of spin-orbit 
misalignment in the XO-3 system (M - 12.5 Mj, P - li.l days, 
e - 0.29) still awaits confirmation as discussed above. 
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The results of the spin-orbit alignment measurements for the 
HD 17156 system can be compared with the prediction of the 
planet scattering models. A range o f alignments can be the out- 
come of planet-planet scattering (iMarzari & Weidenschillin j, 
I2002h . Therefore, our indication for coplanarity does not exclude 
planet-planet scattering in the HD 17156 system. A larger num- 
ber of transiting planets with significant eccentricities have to be 
discovered and characterized to get more conclusive inferences. 

We also searched for stellar companions using adaptive op- 
tics, to test the hypothesis of Kozai mechanism to explain the 
large eccentricity of HD 17156 b. We did not detect compan- 
ions within 1 000 AU, and our detection limits allowed us to 
exclude main sequence companions with projected separations 
from about 150 to 1 000 AU. This result makes unlikely the 
occurrence of a companion inducing Kozai eccentricity oscil- 
lations on the planet, but this possibility can not be yet com- 
pletely rule out (companions at small projected separation and 
faint white dwarfs and brown dwarfs companion still possible). 
Continuation of radial velocity and photometric monitoring will 
allow a more complete view on the existence of additional com- 
panions at small separations. 
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